Comparative Analysis of Driver Mutations and Transcriptomes in Papillary Thyroid Cancer by Region of Residence in South Korea.

BACKGRUOUND: Radiation exposure is a well-known risk factor for papillary thyroid cancer (PTC). South Korea has 24 nuclear reactors in operation; however, no molecular biological analysis has been performed on patients with PTC living near nuclear power plants. METHODS: We retrospectively included patients with PTC (n=512) divided into three groups according to their place of residence at the time of operation: inland areas (n=300), coastal areas far from nuclear power plants (n=134), and nuclear power plant areas (n=78). After propensity score matching (1:1:1) by age, sex, and surgical procedure, the frequency of representative driver mutations and gene expression profiles were compared (n=50 per group). Epithelial-mesenchymal transition (EMT), BRAF, thyroid differentiation, and radiation scores were calculated and compared. RESULTS: No significant difference was observed in clinicopathological characteristics, including radiation exposure history and the frequency of incidentally discovered thyroid cancer, among the three groups. BRAFV600E mutation was most frequently detected in the groups, with no difference among the three groups. Furthermore, gene expression profiles showed no statistically significant difference. EMT and BRAF scores were higher in our cohort than in cohorts from Chernobyl tissue bank and The Cancer Genome Atlas Thyroid Cancer; however, there was no difference according to the place of residence. Radiation scores were highest in the Chernobyl tissue bank but exhibited no difference according to the place of residence. CONCLUSION: Differences in clinicopathological characteristics, frequency of representative driver mutations, and gene expression profiles were not observed according to patients' region of residence in South Korea.

Immune response and mesenchymal transition of papillary thyroid carcinoma reflected in ultrasonography features assessed by radiologists and deep learning.

INTRODUCTION: Ultrasonography (US) features of papillary thyroid cancers (PTCs) are used to select nodules for biopsy due to their association with tumor behavior. However, the molecular biological mechanisms that lead to the characteristic US features of PTCs are largely unknown. OBJECTIVES: This study aimed to investigate the molecular biological mechanisms behind US features assessed by radiologists and three convolutional neural networks (CNN) through transcriptome analysis. METHODS: Transcriptome data from 273 PTC tissue samples were generated and differentially expressed genes (DEGs) were identified according to US feature. Pathway enrichment analyses were also conducted by gene set enrichment analysis (GSEA) and ClusterProfiler according to assessments made by radiologists and three CNNs - CNN1 (ResNet50), CNN2 (ResNet101) and CNN3 (VGG16). Signature gene scores for PTCs were calculated by single-sample GSEA (ssGSEA). RESULTS: Individual suspicious US features consistently suggested an upregulation of genes related to immune response and epithelial-mesenchymal transition (EMT). Likewise, PTCs assessed as positive by radiologists and three CNNs showed the coordinate enrichment of similar gene sets with abundant immune and stromal components. However, PTCs assessed as positive by radiologists had the highest number of DEGs, and those assessed as positive by CNN3 had more diverse DEGs and gene sets compared to CNN1 or CNN2. The percentage of PTCs assessed as positive or negative concordantly by radiologists and three CNNs was 85.6% (231/273) and 7.1% (3/273), respectively. CONCLUSION: US features assessed by radiologists and CNNs revealed molecular biologic features and tumor microenvironment in PTCs.

Impact of thyroid hormone replacement on the risk of second cancer after thyroidectomy: a Korean National Cohort Study.

We aimed to investigate the effect of thyroid hormone administration on the risk of second primary cancer in patients who underwent thyroidectomy for differentiated thyroid cancer. Data were extracted from the medical billing data of the Health Insurance Review and Assessment Service in South Korea. Patients between 19 and 80 years old who underwent thyroid surgery at least once between January 2009 and June 2020 were included. Data of patients with second primary cancer and control patients with matched age, sex, operation date, and follow-up duration were extracted at a ratio of 1:4. A nested case-control analysis was performed to exclude length bias to confirm the correlation between the duration of thyroid hormone administration, dose, and incidence of second primary cancer. Of the 261,598 patients who underwent surgery for thyroid cancer included in the study, 11,790 with second primary cancer and 47,160 without second primary cancer were matched. The average dose of thyroid hormone increased the adjusted odds ratio (OR) for both low (≤ 50 µg, OR 1.29, confidence interval (CI) 1.12-1.48) and high (< 100 µg, OR 1.24, CI 1.12-1.37) doses. Analyzing over time, the adjusted OR of second primary cancer increased, especially in short (≤ 1 year) (OR 1.19; CI 1.06-1.34) and long (> 5 years) duration (OR 1.25; CI 1.10-1.41). In conclusion, insufficient and excessive thyroid hormone replacement might be linked to increased second primary cancer in patients who underwent thyroidectomy for differentiated thyroid cancer.

Impact of thyroid cancer on the cancer risk in patients with non-alcoholic fatty liver disease or dyslipidemia.

The raised prevalence of obesity has increased the incidence of obesity-related metabolic diseases such as dyslipidemia (DL) and non-alcoholic fatty liver disease (NAFLD), along with the development and progression of various types of cancer, including thyroid cancer. In this study, we investigated whether thyroid cancer in patients with DL and NAFLD could be a risk factor for other cancers. To achieve our goal, we generated two independent cohorts from our institution and from the National Health Insurance System in South Korea. Based on the ICD-10 code, we conducted exact matching (1:5 matching) and estimated the overall risk of thyroid cancer for other cancers in patients with DL or NAFLD. Univariate and multivariate analyses showed that the hazard ratio (HR) of thyroid cancer was 2.007 (95% Confidence Interval [CI], 1.597-2.522) and 2.092 (95% CI, 1.546-2.829), respectively in the institutional cohort and 1.329 (95% CI, 1.153-1.533) and 1.301 (95% CI, 1.115-1.517), respectively in the nationwide cohort. Risk analysis revealed a significant increase in the HR in lip, tongue, mouth, lung, bone, joint, soft tissue, skin, brain, male cancers and lymphoma after thyroid cancer occurred. Thyroid cancer in patients with DL or NAFLD might be a valuable factor for predicting the development of other cancers.

Heteroepitaxial Growth of B5 -Site-Rich Ru Nanoparticles Guided by Hexagonal Boron Nitride for Low-Temperature Ammonia Dehydrogenation.

Ruthenium is one of the most active catalysts for ammonia dehydrogenation and is essential for the use of ammonia as a hydrogen storage material. The B5 -type site on the surface of ruthenium is expected to exhibit the highest catalytic activity for ammonia dehydrogenation, but the number of these sites is typically low. Here, a B5 -site-rich ruthenium catalyst is synthesized by exploiting the crystal symmetry of a hexagonal boron nitride support. In the prepared ruthenium catalyst, ruthenium nanoparticles are formed epitaxially on hexagonal boron nitride sheets with hexagonal planar morphologies, in which the B5 sites predominate along the nanoparticle edges. By activating the catalyst under the reaction condition, the population of B5 sites further increases as the facets of the ruthenium nanoparticles develop. The electron density of the Ru nanoparticles also increases during catalyst activation. The synthesized catalyst shows superior catalytic activity for ammonia dehydrogenation compared to previously reported catalysts. This work demonstrates that morphology control of a catalyst via support-driven heteroepitaxy can be exploited for synthesizing highly active heterogeneous catalysts with tailored atomic structures.

Non-genomic activation of the AKT-mTOR pathway by the mitochondrial stress response in thyroid cancer.

Cancer progression is associated with metabolic reprogramming and causes significant intracellular stress; however, the mechanisms that link cellular stress and growth signalling are not fully understood. Here, we identified a mechanism that couples the mitochondrial stress response (MSR) with tumour progression. We demonstrated that the MSR is activated in a significant proportion of human thyroid cancers via the upregulation of heat shock protein D family members and the mitokine, growth differentiation factor 15. Our study also revealed that MSR triggered AKT/S6K signalling by activating mTORC2 via activating transcription factor 4/sestrin 2 activation whilst promoting leucine transporter and nutrient-induced mTORC1 activation. Importantly, we found that an increase in mtDNA played an essential role in MSR-induced mTOR activation and that crosstalk between MYC and MSR potentiated mTOR activation. Together, these findings suggest that the MSR could be a predictive marker for aggressive human thyroid cancer as well as a useful therapeutic target.

Cell non-autonomous effect of hepatic growth differentiation factor 15 on the thyroid gland.

The thyroid gland plays an essential role in the regulation of body energy expenditure to maintain metabolic homeostasis. However, to date, there are no studies investigating the morphological and functional changes of the thyroid gland due to mitochondrial stress in metabolic organs such as the liver. We used data from the Genotype-Tissue Expression portal to investigate RNA expression patterns of the thyroid gland according to the expression of growth differentiation factor 15 (GDF15) such as the muscles and liver. To verify the effect of hepatic GDF15 on the thyroid gland, we compared the morphological findings of the thyroid gland from liver-specific GDF15 transgenic mice to that of wild type mice. High GDF15 expression in the muscles and liver was associated with the upregulation of genes related to hypoxia, inflammation (TGF-α via NFκB), apoptosis, and p53 pathway in thyroid glands. In addition, high hepatic GDF15 was related to epithelial mesenchymal transition and mTORC1 signaling. Electron microscopy for liver-specific GDF15 transgenic mice revealed short mitochondrial cristae length and small mitochondrial area, indicating reduced mitochondrial function. However, serum thyroid stimulating hormone (TSH) level was not significantly different. In our human cohort, those with a high serum GDF15 level showed high fasting glucose, alanine transaminase, and alkaline phosphatase but no difference in TSH, similar to the data from our mice model. Additionally, high serum GDF15 increased the risk of lymph node metastasis to lateral neck. The hepatic GDF15 affected thyroid morphogenesis via a TSH-independent mechanism, affecting aggressive features of thyroid cancers.

Long Non-Coding RNA-Based Functional Prediction Reveals Novel Targets in Notch-Upregulated Ovarian Cancer.

Notch signaling is a druggable target in high-grade serous ovarian cancers; however, its complexity is not clearly understood. Recent revelations of the biological roles of lncRNAs have led to an increased interest in the oncogenic action of lncRNAs in various cancers. In this study, we performed in silico analyses using The Cancer Genome Atlas data to discover novel Notch-related lncRNAs and validated our transcriptome data via NOTCH1/3 silencing in serous ovarian cancer cells. The expression of novel Notch-related lncRNAs was down-regulated by a Notch inhibitor and was upregulated in high-grade serous ovarian cancers, compared to benign or borderline ovarian tumors. Functionally, Notch-related lncRNAs were tightly linked to Notch-related changes in diverse gene expressions. Notably, genes related to DNA repair and spermatogenesis showed specific correlations with Notch-related lncRNAs. Master transcription factors, including EGR1, CTCF, GABPα, and E2F4 might orchestrate the upregulation of Notch-related lncRNAs, along with the associated genes. The discovery of Notch-related lncRNAs significantly contributes to our understanding of the complex crosstalk of Notch signaling with other oncogenic pathways at the transcriptional level.

Simultaneous Expression of Long Non-Coding RNA FAL1 and Extracellular Matrix Protein 1 Defines Tumour Behaviour in Young Patients with Papillary Thyroid Cancer.

We investigated the regulatory mechanism of FAL1 and unravelled the molecular biological features of FAL1 upregulation in papillary thyroid cancer (PTC). Correlation analyses of FAL1 and neighbouring genes adjacent to chromosome 1q21.3 were performed. Focal amplification was performed using data from copy number alterations in The Cancer Genome Atlas (TCGA) database. To identify putative transcriptional factors, PROMO and the Encyclopaedia of DNA Elements (ENCODE) were used. To validate c-JUN and JUND as master transcription factors for FAL1 and ECM1, gene set enrichment analysis was performed according to FAL1 and ECM1 expression. Statistical analyses of the molecular biological features of FAL1- and ECM1-upregulated PTCs were conducted. FAL1 expression significantly correlated with that of neighbouring genes. Focal amplification of chromosome 1q21.3 was observed in ovarian cancer but not in thyroid carcinoma. However, PROMO suggested 53 transcription factors as putative common transcriptional factors for FAL1 and ECM1 simultaneously. Among them, we selected c-JUN and JUND as the best candidates based on ENCODE results. The expression of target genes of JUND simultaneously increased in FAL1- and ECM1-upregulated PTCs, especially in young patients. The molecular biological features represented RAS-driven PTC and simultaneously enriched immune-related gene sets. FAL1 and ECM1 expression frequently increased simultaneously and could be operated by JUND. The simultaneous upregulation might be a potential diagnostic and therapeutic target for RAS-driven PTC.

Impact of Dyslipidemia on the Risk of Second Cancer in Thyroid Cancer Patients: A Korean National Cohort Study.

BACKGROUND: Studies have shown that radioactive iodine therapy (RAIT) affects the development of second cancer in thyroid cancer patients. The impact of other factors, such as dyslipidemia are not clear. METHODS: A retrospective analysis of thyroid cancer patients with a 1,251,913 person-year follow-up was conducted using data from the Health Insurance Review and Assessment database in South Korea from January 2008 to December 2018. We investigated factors related to second cancer development using a nested case-control analysis to avoid length bias. RESULTS: The overall risk of developing second cancer was higher in thyroid cancer patients than in the general population [standardized incidence ratio, 3.34; 95% confidence interval (CI) 3.30-3.39]. Second cancer incidence was higher in patients who received RAIT than in those who did not [odds ratio (OR) 1.130; 95% CI 1.094-1.169]. Moreover, the risk of second cancer was higher in patients with dyslipidemia than in those without dyslipidemia (OR 1.265; 95% CI 1.223-1.309). After adjustment for RAIT, the incidence of a second cancer was higher in patients with dyslipidemia than in those without dyslipidemia (OR 1.262; 95% CI 1.221-1.306). CONCLUSIONS: The risk of second cancer development in patients with thyroid cancer appears to be high. Dyslipidemia may be associated with an increased risk of several types of second cancers.

Top-Down Syntheses of Nickel-Based Structured Catalysts for Hydrogen Production from Ammonia.

We report the fabrication and catalytic performance evaluation of highly active and stable nickel (Ni)-based structured catalysts for ammonia dehydrogenation with nearly complete conversion using nonprecious metal catalysts. Low-temperature chemical alloying (LTCA) followed by selective aluminum (Al) dealloying was utilized to synthesize foam-type structured catalysts ready for implementation in commercial-scale catalytic reactors. The crystalline phases of Ni-Al alloy (NiAl3, Ni2Al3, or both) in the near-surface layer were controlled by tuning the alloying time. The best-performing catalyst was obtained from a Ni foam substrate with a NiAl3/Ni2Al3 overlayer synthesized by LTCA at 400 °C for 20 h. The developed Ni catalyst exhibited an activity enhancement of 10-fold over the nontreated Ni foam and showed outstanding activities of 15 800 molH2molNi-1h-1 (TOF: 4.39 s-1) and 19 978 molH2molNi-1h-1 (TOF: 5.55 s-1) at 550 and 600 °C, respectively. This performance is unprecedented compared with previously reported Ni-based ammonia cracking catalysts with higher-end performance (TOFs of 0.08-1.45 s-1 at 550 °C). Moreover, this catalyst showed excellent stability for 100 h at 600 °C while discharging an extremely low NH3 concentration of 1034 ppm. The NH3 concentration in the exhaust gas was further reduced to 690 and 271 ppm at 700 and 800 °C, respectively, while no deactivation was observed at these elevated temperatures. Through material characterizations, we clarified that controlling the degree of Al alloying in the outermost layer of Ni is a crucial factor in determining the activity and stability because residual Al possibly modifies the electronic structure of Ni for enhanced activity as well as transforming to acidic alumina for increased intrinsic activity and stability.

Artificial intelligence to predict the BRAFV600E mutation in patients with thyroid cancer.

PURPOSE: To investigate whether a computer-aided diagnosis (CAD) program developed using the deep learning convolutional neural network (CNN) on neck US images can predict the BRAFV600E mutation in thyroid cancer. METHODS: 469 thyroid cancers in 469 patients were included in this retrospective study. A CAD program recently developed using the deep CNN provided risks of malignancy (0-100%) as well as binary results (cancer or not). Using the CAD program, we calculated the risk of malignancy based on a US image of each thyroid nodule (CAD value). Univariate and multivariate logistic regression analyses were performed including patient demographics, the American College of Radiology (ACR) Thyroid Imaging, Reporting and Data System (TIRADS) categories and risks of malignancy calculated through CAD to identify independent predictive factors for the BRAFV600E mutation in thyroid cancer. The predictive power of the CAD value and final multivariable model for the BRAFV600E mutation in thyroid cancer were measured using the area under the receiver operating characteristic (ROC) curves. RESULTS: In this study, 380 (81%) patients were positive and 89 (19%) patients were negative for the BRAFV600E mutation. On multivariate analysis, older age (OR = 1.025, p = 0.018), smaller size (OR = 0.963, p = 0.006), and higher CAD value (OR = 1.016, p = 0.004) were significantly associated with the BRAFV600E mutation. The CAD value yielded an AUC of 0.646 (95% CI: 0.576, 0.716) for predicting the BRAFV600E mutation, while the multivariable model yielded an AUC of 0.706 (95% CI: 0.576, 0.716). The multivariable model showed significantly better performance than the CAD value alone (p = 0.004). CONCLUSION: Deep learning-based CAD for thyroid US can help us predict the BRAFV600E mutation in thyroid cancer. More multi-center studies with more cases are needed to further validate our study results.

Liver X Receptor ß Related to Tumor Progression and Ribosome Gene Expression in Papillary Thyroid Cancer.

BACKGROUND: Intracellular lipid deposition has been reported in thyroid glands in obese animal and human. To understand the regulatory mechanism of lipid metabolism in thyroid cancer, we investigated the expression status of liver X receptor (LXR) and analyzed its clinicopathological characteristics and molecular biological features. METHODS: Expression status of LXR and its transcriptional targets in human cancers were analyzed using The Cancer Genome Atlas (TCGA). The gene-sets related to high LXRß expression was investigated by gene set enrichment analysis (GSEA) using Kyoto Encyclopedia of Genes and Genomes (KEGG) signaling pathways and gene ontology biologic process. Quantitative reverse transcription polymerase chain reaction was performed in thyroid cancer samples using our validation cohort. RESULTS: In contrast to low expression of LXRα, LXRß was highly expressed in thyroid cancer compared to the other types of human cancers. High LXRß expression was correlated with the expression of LXRß transcriptional targets genes, such as apolipoprotein C1 (APOC1), APOC2, apolipoprotein E (APOE), ATP binding cassette subfamily G member 8 (ABCG8), sterol regulatory elementbinding protein 1c (SREBP1c), and SPOT14. Furthermore, High LXRß expression group indicated poor clinicopathological characteristics and aggressive molecular biological features independently from the drive mutation status. Mechanistically, high LXRß expression was coordinately related to ribosome-related gene sets. CONCLUSION: The mechanistic link between LXRß and ribosomal activity will be addressed to develop new diagnostic and therapeutic targets in thyroid cancers.

Detailed characterization of metastatic lymph nodes improves the prediction accuracy of currently used risk stratification systems in N1 stage papillary thyroid cancer.

OBJECTIVE: The characteristics of metastatic lymph nodes (MLNs) have been investigated as important predictors of recurrence and progression in papillary thyroid cancer (PTC). However, clinically applicable risk stratification systems are limited to the assessment of size and number of MLNs. This study investigated the predictive value of detailed characteristics of MLNs in combination with currently used risk stratification systems. DESIGN AND METHODS: We retrospectively characterized 2811 MLNs from 9014 harvested LNs of 286 patients with N1 PTC according to the maximum diameter of MLN (MDLN), maximum diameter of metastatic focus (MDMF), ratio of both diameters (MDMFR), lymph node ratio (LNR, number of MLNs/number of total harvested LNs), presence of extranodal extension (ENE), desmoplastic reaction (DR), cystic component, and psammoma body. RESULTS: Factors related to the size and number of MLNs were associated with increased risk of recurrence and progression. Extensive presence of ENE (>40%) and DR (≥50%) increased the risk of recurrence/progression. The combination of MDLN, LNR, ENE, and DR had the highest predictive value among MLN characteristics. Combination of these parameters with ATA risk stratification or 1-year response to therapy improved the predictive power for recurrence/progression from a Harrell's C-index of 0.781 to 0.936 and 0.867 to 0.960, respectively. CONCLUSIONS: The combination of currently used risk stratification systems with detailed characterization of MLNs may improve the predictive accuracy for recurrence/progression in N1 PTC patients.

Effects of Oxytocin on Cell Proliferation in a Corticotroph Adenoma Cell Line.

BACKGROUND: Oxytocin (OXT) has been reported to act as a growth regulator in various tumor cells. However, there is a paucity of data on the influence of OXT on cell proliferation of corticotroph adenomas. This study aimed to examine whether OXT affects cell growth in pituitary tumor cell lines (AtT20 and GH3 cells) with a focus on corticotroph adenoma cells. METHODS: Reverse transcription polymerase chain reaction and enzyme-linked immunosorbent assay were conducted with AtT20 cells to confirm the effects of OXT on hormonal activity; flow cytometry was used to assess changes in the cell cycle after OXT treatment. Moreover, the impact of OXT on proliferating cell nuclear antigen (PCNA), nuclear factor κB, and mitogen-activated protein kinase signaling pathway was analyzed by Western blot. RESULTS: OXT treatment of 50 nM changed the gene expression of OXT receptor and pro-opiomelanocortin within a short time. In addition, OXT significantly reduced adrenocorticotropic hormone secretion within 1 hour. S and G2/M populations of AtT20 cells treated with OXT for 24 hours were significantly decreased compared to the control. Furthermore, OXT treatment decreased the protein levels of PCNA and phosphorylated extracellular-signal-regulated kinase (P-ERK) in AtT20 cells. CONCLUSION: Although the cytotoxic effect of OXT in AtT20 cells was not definite, OXT may blunt cell proliferation of corticotroph adenomas by altering the cell cycle or reducing PCNA and P-ERK levels. Further research is required to investigate the role of OXT as a potential therapeutic target in corticotroph adenomas.

Clinical Value of Lymph Node Ratio Integration with the 8th Edition of the UICC TNM Classification and 2015 ATA Risk Stratification Systems for Recurrence Prediction in Papillary Thyroid Cancer.

Recently, the 2015 American Thyroid Association (ATA) risk stratification and the 8th edition of the American Joint Committee on Cancer/Union for International Cancer Control (AJCC/UICC) TNM staging system were released. This study was conducted to assess the clinical value of the lymph node ratio (LNR) as a predictor of recurrence when integrated with these newly released stratification systems, and to compare the predictive accuracy of the modified systems with that of the newly released systems. The optimal LNR threshold value for predicting papillary thyroid cancer (PTC) recurrence was 0.17857 using the Contal and O'Quigley method. The 8th edition of the AJCC/UICC TNM staging system with the LNR and the 2015 ATA risk stratification system with the LNR were significant predictors of recurrence. Furthermore, calculation of the proportion of variance explained (PVE), the Akaike information criterion (AIC), Harrell's c index, and the incremental area under the curve (iAUC) revealed that the 8th edition of the TNM staging system with the LNR, and the 2015 ATA risk stratification system with the LNR, showed the best predictive performance. Integration of the LNR with the TNM staging and the ATA risk stratification systems should improve prediction of recurrence in patients with PTC.

Peripheral location and infiltrative margin predict invasive features of papillary thyroid microcarcinoma.

OBJECTIVE: Tumor location in papillary thyroid microcarcinoma (PTMC) might determine tumor outgrowth from the thyroid gland. However, the clinical implications of tumor location and minimal extrathyroid extension (mETE) have not been well elucidated. We aimed to investigate the relationship between tumor location and mETE to predict the aggressiveness of PTMC. METHODS: A total of 858 patients with PTMC were grouped according to tumor location on ultrasonography: central (cPTMC) and peripheral PTMC (pPTMC). PTMC without mETE (PTMC-mETE(-)) was divided further according to margin shape: encapsulated (E-) or infiltrative (I-). To understand the molecular biologic characteristics of PTMC presenting with an I-margin and mETE, transcriptome data from TCGA-THCA were analyzed using Gene Set Enrichment Analysis (GSEA). RESULTS: pPTMC (n = 807, 94.1%) accounted for the majority of cases; mETE was identified only in pPTMC (403/807; 49.9%). pPTMC-mETE(+) showed aggressive clinical characteristics that increased the odds ratio (OR) for lymph node metastasis (LNM). Interestingly, subgroup analysis of PTMC-mETE(-) revealed that the I-margin also increased the OR for LNM, independent of other clinical factors. GSEA of TCGA-THCA data suggested coordinated upregulation of genes related to epithelial-mesenchymal transition (EMT) in PTC with mETE. Immunohistochemical staining for laminin subunit gamma 2 (LAMC2), CD59, E-cadherin and vimentin showed that these markers of EMT were associated with progressive changes in E-margin PTMC-mETE(-), I-margin PTMC-mETE(-) and pPTMC-mETE(+). CONCLUSION: mETE related to peripheral location of PTMC is an important predictor of tumor invasiveness, as is the I-margin, which presents with EMT features similar to mETE. I-margin PTMC-mETE(-) and pPTMC-mETE(+) might reflect the pattern of invasive PTMC.

Whole Exome Sequencing Identifies a Novel Hedgehog-Interacting Protein G516R Mutation in Locally Advanced Papillary Thyroid Cancer.

Locally advanced thyroid cancer exhibits aggressive clinical features requiring extensive neck dissection. Therefore, it is important to identify changes in the tumor biology before local progression. Here, whole exome sequencing (WES) using tissues from locally advanced papillary thyroid cancer (PTC) presented a large number of single nucleotide variants (SNVs) in the metastatic lymph node (MLN), but not in normal tissues and primary tumors. Among those MLN-specific SNVs, a novel HHIP G516R (G1546A) mutation was also observed. Interestingly, in-depth analysis for exome sequencing data from the primary tumor presented altered nucleotide 'A' at a very low frequency indicating intra-tumor heterogeneity between the primary tumor and MLN. Computational prediction models such as PROVEAN and Polyphen suggested that HHIP G516R might affect protein function and stability. In vitro, HHIP G516R increased cell proliferation and promoted cell migration in thyroid cancer cells. HHIP G516R, a missense mutation, could be a representative example for the intra-tumor heterogeneity of locally advanced thyroid cancer, which can be a potential future therapeutic target for this disease.